Sigma models on curved space and bubble refraction in doped 

antiferromagnets 

J.M. Speight 

School of Mathematics, University of Leeds, Leeds LS2 9JT, UK 

Abstract 

The dynamics of bubble solitons in two-dimensional isotropic antiferromagnets, inhomogeneously doped 
so that the exchange integral J becomes position dependent, is studied. In the usual continuum approxi- 
mation, the system reduces to a nonlinear sigma model on a spacetime whose geometry depends on J{x). 
It is shown, both within the geodesic approximation, and by appealing to field theoretic conservation laws, 
that a bubble incident on a domain-wall inhomogeneity undergoes refraction in accordance with Snell's law, 
J -1 being identified with the refractive index, and that sufficiently oblique impacts result in total internal 
reflexion. Possible applications of this phenomenon to the construction of bubble lenses and bubble guides 
(in analogy with fibre-optic cables) are considered. 



1 Introduction 

Topological solitons have important applications in many branches of physics, especially condensed matter 
and high energy physics. It often happens that the same soliton can be given completely different physical 
interpretations in these two disciplines, which has led in recent years to a very fruitful exchange of ideas and 
techniques. One of the most surprising and least exploited of these soliton correspondences is that between the 
isotropic Heisenberg antiferromagnet and the relativistic 0(3) sigma model. In this paper we will exploit this 
correspondence to show, using techniques from particle theory and differential geometry, that the trajectories 
of bubble-like solitons in inhomogeneously doped two-dimensional antiferromagnets experience refraction in 
exact analogy with Snell's law of refraction in geometric optics. 

The isotropic Heisenberg antiferromagnet consists of an infinite square lattice of classical spin variables 
Sij G R 3 , with \Sij\ = s 2 , evolving according to the law 

dS i3 „ dH 



S « x ^rT' H :=£^[2s 2 + S ij -(S i . j+1 + S i+1J -)] (1) 
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where r is time and J > is a constant called the exchange integral. Note that since J > 0, H is minimized 
when each spin anti-aligns with its nearest neighbours. It is a curious fact that this dynamical system, first 
order in time and with an obvious preferred reference frame, can be described, in the continuum limit, by a 
Lorentz invariant PDE which is second order in time, namely 

nx^n = 0, ^ — -J^— + — J (2) 

where n : R 2+1 — ► R 3 has |n| 2 = 1, and t is a rescaled time variable. Note that J5J) holds if and only if is 
parallel to n, and hence 

^n-(n-^n)n = 0, (3) 

which is the more familar form of the 0(3) sigma model (or wave map) equation. The connexion between n 
and Sy is subtle, and will be explained in the next section. Roughly speaking, one should imagine the field 
n sampled on a square lattice, partitioned into black and white sublattices, chess-board fashion. The white 
spins evolve as sn(t), the black as —sn(t). 



Imagine now that the antiferromagnet is inhomogcneously doped, so that the exchange integral becomes a 
function of position, J(x,y). We will argue that the spin dynamics is now described by a field n obeying 



2„ /P2„ tf,, , 

0, (4) 



d 2 n 2 / d 2 n d 2 n 



which is the 0(3) sigma model equation on spacetime R x M 2 equipped with an inhomogeneous Lorentzian 
metric 

V = dt 2 - -j^-^idx'+dy 2 ). (5) 

Note that the static field equation is independent of J, so the usual Belavin-Polyakov lumps of the homogeneous 
system 0, which we shall call bubbles in this context (in analogy with magnetic bubbles in ferromagnets), 
carry over unchanged to the doped system. In this paper we make a detailed study of the dynamics of a single 
bubble in various inhomogeneous antiferromagnets, primarily within the geodesic approximation of Manton 
0]. The main calculational task in this approach is to compute the metric on the one-bubble moduli space. 
This metric does depend on J(x,y), and we obtain exact expressions for it for several interesting choices of 
this function. 

Our main conclusion is that bubble trajectories experience refraction, in exact analogy with geometric 
optics, J{x,y)^ 1 playing the role of the refractive index of the medium. The incident and exit angles of a 
bubble crossing a domain-wall J-inhomogeneity, for example, are related by Snell's law, and total internal 
reflexion occurs if the impact is sufficiently oblique. This key observation is rederived without reference to the 
geodesic approximation by appealing to field theoretic conservation laws. We go on to study simple prototype 
bubble "lenses" and bubble guides, which use total internal reflexion to guide bubble trajectories just as fibre 
optic cables guide light beams. 

The rest of the paper is structured as follows. In section [5] we derive the PDE (@J from 0J. In section |21 
the geodesic approximation to the dynamics of a single bubble is set up. The interaction of a bubble with a 
domain wall is studied inland Snell's law of refraction obtained. Sections |5] and [5] consider simple lenses and 
bubble guides respectively, while presents some concluding remarks. 



2 Doped antiferromagnets in the continuum limit 

The spin lattice of interest is defined by (JJJ, but with J replaced by a position dependent exchange integral, 

J ij , SO 

= ./, ; S, x (S l+1 ,, + S, , ; + S id+1 + S, ; ,:. (6) 

Following pQ, to obtain a well-defined continuum limit for this system we, think of it as a lattice of dimers, that 
is, spin pairs. First, partition the square lattice chess-board fashion into black and white sublattices, with site 
white (black) if i — j is even (odd). Associate to each white site (i, j) the dimer consisting of (Sy, Sj.j+i). 
The white sites can be labelled by a pair of (unconstrained) integers (a,/3) = [{i + j)/2, (i — j)/2). We then 
call the two spins of the dimer A a( g = Sij and B a £j = Sjj+i. Note that, while A lives on the white lattice 
and B on the black, both sets of variables are labelled with reference to the white lattice. Their evolution is 
governed by 



dA 



a0 



dr 
dr 



■J<x-/3,a+p(Ba,0-l + B a /3 + B Q _i j( g + B a _i i( g_i) (7) 

~Jot—8.a +/3+i( a q+i,/3 + A Q+ i !/3+ i + A a> p + i + A Qj/3 ). (8) 



We now think of the original lattice as having spacing e > 0, small, so that the (a, (3) lattice has 

spacing 5 — \[2e, and assume that A a ^, B Qi a and Jij have well defined continuum limits A(?7,£), B(?7,^) 
and J(r],t;), where ij — aS, £ = (36, so that A a +i )( g = A(rj,£) + (JA,, -I- ^5 2 A V11 + • • • , and so on. Note that 

JfaiO = Ja-P.a+P, SO 

J a -p, a +0+i = J(»7+^,£+~<5) = J(^g) + £(J, + J c ) + i.(J IW + J £C + 2J rrf ) + .-. . (9) 
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Substituting into Q, JHJ and discarding terms of order S 3 or higher, one obtains 

Be) +S 2 (B, m + B e 



2sSA t 
2s6B t 



-JA x [4B - 26{B V 
-JB x [4A + 2<5(A ?? 



A £ ) 



5 (B 7p] 
S (A^ 



5 ?? 

A 55 



■ B, £ )] 



-B x [2£(J„ + J C )A + <5 2 (J, ; + J e )(A„ + A c ) + —{i m 



J« + 2J^)A], 



where t — 2sSt is a rescaled time variable. 

These equations simplify somewhat if we introduce the auxilliary fields 



1 

2s 



(A + B) 



1 

2s 



(A-B). 



(10) 



(11) 



(12) 



Note that m • n = and |m| 2 + |n| 2 = 1. Since we expect neighbouring spins to almost anti-align, m should 
be small, so we will assume, consistently with (ll()|) . i[ll|l. that |m| = 0(5), whence |n| = 1 + 0(S 2 ). Then, 
again neglecting terms of order higher than <5 2 , one finds 



-(d v + <%)[J m x n] + -[2J n x (n^ 



Sn t = 4J m x n - if J n x (n,, 
Equation Ijl4(l may be solved for m, 



•n*)-*(J„ + J*)mxn-— (J, 



J 4 )n x (n^ - 
- Jj)n x (n ? , 



•"€)] 
-n 6 ). 



8 A 
- -n 
4 L J 



x n 4 



0(^ 2 



(13) 
(14) 

(15) 



which may then be substituted into equation (|13|l . One finds that all the terms involving derivatives of J 
miraculously cancel, and 

nxn lf = J 2 nx(n„ + n K )+0(i). (16) 

To leading order, therefore, n is governed by the 0(3) sigma model equation on 1 x I 2 with metric dt 2 — 
(drf + d£ 2 )/J(7/, £) 2 , which is Q up to relabelling. It is slightly more natural, however, to define continuum 
position variables x — ie, y — je aligned with the original spin lattice, rather than r), £ aligned with the white 
sublattice. Since the x,y coordinates are obtained from r], £ by a rotation (through 45°), d 2 + d 2 = d 2 + d 2 , 
so we again obtain but with J(x,y) = Jij. To reconstruct the actual spin dynamics, given a solution of 
(J3J, we must choose e > small, construct the auxilliary field m from (|15fl . 



m=^[inXn t -„ 2/ ], 



(17) 



extract A and B, then sample these on the white and black sublattices, respectively, of the lattice of spacing e. 

The above derivation is closely modelled on that of Komineas and Papanicolaou pQ, so we have omitted 
some of the details of the calculation. The main difference is that we have included the extra terms due to the 
position dependence of J. In the end, all these extra terms cancel, to leading order, and we are left with the 
continuum PDE one would have first guessed on naive grounds. This is somewhat surprising, and it would be 
satisfying to have a more direct argument for it than the rather complicated derivation given here. 



3 The geodesic approximation 

Equation {2J is the variational equation for the action 

S=^Jdtdxdyy/\^\d lJl n-d v nr]> iV = J dt (T — V) 



where 



T = \S dxdy l^W 



On 



dt 



V = - / dxdy 



9n 



Ox 



On 



dy 



(18) 



(19) 
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are identifed as the kinetic and potential energy functionals respectively. Note that while T depends on J, V 
does not. Hence Belavin's and Polyakov's analysis of the static homogeneous model carries over unchanged 
to the doped system. There is a topological lower energy bound 

V[n] > 4ttH, (20) 

m being the topological degree of the map n : M 2 U {00} — » S 2 . Let us define complex coordinates z = x + iy 
on the spatial plane and u = (n\ + ma)/(l — TI3) on the target sphere (the latter being the image of n under 
stereographic projection from (0,0, 1) to the equatorial plane). Then the bound (|2()|> for m > is attained 
if and only if u(z) is a rational map of algebraic degree m. Since such rational maps minimize V (globally 
within their homotopy class, in fact) they are automatically stable static solutions of the model. If we choose 
the boundary value of n at spatial infinity to be (0, 0, 1), the general degree 1 bubble is 

u ( z ) = X^ 1 e~ l,p (z - a) =: u a (z; x, ip, a), (21) 

X G (0, 00), ip £ [0, 27r] and a = a\ + ia% G C being constants interpreted as the bubble's width, internal phase 
and position respectively. 

Following Ward [H] and Leese 0, we will study the low energy dynamics of a single bubble within the 
geodesic approximation of Manton [3]. The idea is that Mi, the four dimensional space of static bubbles, 
parametrized by x, ip, a\, 02, which we will denote collectively as q l , is a flat valley bottom in the space of all 
degree 1 maps R 2 — > S 2 , on which V attains its minimum value of Air. If a bubble is given a small amount of 
kinetic energy, its subsequent motion is confined close to Mi by conservation of energy E — T + V , suggesting 
that a collective coordinate approximation wherein the motion is constrained to M 1 for all time is sensible. So 
we assume that 

u(z,t) = x(t)e-**<-*\z-a(t)), (22) 

substitute Q22JI into S, and obtain variational equations for x,4>i a - Note that V = 47r, constant, for all fields 
of the form (|22|l . and T is quadratic in time derivatives, so the action reduces to 

c 1 f -i-j r\ [ dxdy A du du , oq x 

2 J J J{x,y) 2 (1 + \u (z;q)\ 2 ) 2 dq 1 dqi 

which is the action for geodesic motion on Mi with respect to the metric g = gij(q)dq l dq J . The conceptual 
framework and validity of this approximation are discussed at length in [7] . 

The metric g, and hence the dynamics, depends strongly on J(x,y). If J is constant, we are studying 
the standard sigma model, and it is found that x an d ip are frozen by infinite inertia (g X x = 94>4> = 00 ) an( i 
g = \-K,J~ 2 da da, so bubbles just travel at constant velocity 0. At the other extreme, if J(x, y) = 1 + x 2 + y 2 , 
we effectively have the sigma model on a round two-sphere, and the bubble dynamics is much richer 8 . This 
choice has J unbounded, which is presumably unphysical, however. If we assume that J remains bounded, 
then x and ip are frozen by an essentially identical argument to Ward's. Without loss of generality, we can set 
ip = 0. The width x remains a free, but frozen, parameter. The metric on Mj, the width Xi phase leaf of 
Mi, is then 

g = f x (a)(da 2 +da 2 ), (24) 



, dzdz Ax 2 f Adudu 1 

fxW = / m I .-si. T^W = / T, I "777— 7^v>> ( 25 ) 



where 

^ J(z,z) 2 (I + x~ 2 \z - a\ 2 ) 2 ~ J (1 + M 2 ) 2 J( X u + a) 2 ' 

the integral of J(xu + a)~ 2 over the standard unit sphere, u being a stereographic coordinate. If J is sufficiently 
simple, f x {a) can be computed explicitly. Note that 

47T 

lim n 9 = -TTyi {da\ + da 2 2 ), (26) 

which is the bubble rest mass, M = 47r, times the inhomogeneous metric on the physical plane. In the limit of 
vanishing width, therefore, bubble trajectories tend to geodesies, not just in Mj, but also in the physical plane. 
One may regard g as a smeared out version, due to the bubble's finite core size, of the metric on physical 
space. 
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4 Domain walls 



The simplest and most experimentally feasible type of inhomogeneous doping is the domain wall. Let us assume 
that J depends only on x, is constant outside a small neighbourhood of x — and rises monotonically from J_ 
to J + as x traverses this neighbourhood. Such a J(x) would arise from doping one end of an antiferromagnet 
but not the other. If the domain wall is very narrow, we may idealize it by a step function, 



J(x) 



x > 
x < 0. 



Then from 125|) we see that f x depends only on ai and 



fxM = -FT 



A- 47r - A- 

J 2 



7 2 



(27) 



(28) 



where A- is the area of 5_ C S 2 , the region on which Ke(\u + a) < 0. The stereographic image of S- is the 
half-plane to the left of the straight line u\ = —ai/x- The corresponding boundary curve dS- C S 2 is a circle 
whose diameter (in the metric space sense) is the length of any great circular arc lying in S- , intersecting dS- 
orthogonally. One such arc has stereographic image {t + iO E C : t G (-co, — cii/x]}- The geometry of the 
situation is depicted in figure ^ Hence 



diam(S'_) 



ai/x 



dt- 



l + t 2 



7r — 2 tan 



Ol 

X ' 



(29) 



To compute the area of any disk in S 2 with diameter diam, we can use spherical polar coordinates with the 
"North pole" placed at the centre of the disk. Then 



A(diam) 

Together, P ^l and pi f yield 

/x(«) = 



2tt />Jj>diam 

dip 
o Jo 



d6 sin 9 — 2n(l ~ cos -diam). 



2tt 
J 2 



Ol 



Vx 2 



2tt 
Ji 



Ol 



Vx 2 



(30) 



(31) 



Note that, despite the discontinuity in J, f x is smooth, and that f x is a decreasing function of a\ with 
lim^—jioo f x (ai) — AttJ± 2 . If J is smoothed around x = 0, these properties, on which we will base all our 
dynamical arguments, persist, although the exact formula (|31|l is lost. 
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To analyze geodesic flow in M*, we note that the Lagrangian governing a(t) is L = hfx( a i){^i + ^1)- The 
time evolution conserves both energy e = L and the momentum conjugate to 02, P2 — f x ( a i)o-2- Since geodesic 
trajectories are independent of initial speed, we may set E = 1. Then 

for all time. Consider a bubble fired obliquely at the domain wall from the left, that is, with ai(0) < 0, 
di(0) > 0, 02(0) 7^ 0. Since / x (ai) is decreasing, a\(t) must have a turning point if \p%\ (hence |a2 (0) |) is too 
large, else (|32|) is violated. Hence a bubble which hits the wall with two large an incidence angle is reflected 
from it, just like a light ray undergoing total internal reflection. In fact, the analogy with refraction is much 
deeper. Let us assume that the incident bubble has velocity t>_(cos#_, sin#_) as t — > —00, — it/2 < 6- < tt/2, 
so that it impinges on the wall, and is transmitted with velocity 1*4. (cos # + , sin# + ) as t — * 00. Then by 
conservation of e and p2, 



2n o 2TT o . 47T . „ 47T 



2 V + — ~~2 V - an d — 2 V + sin#+ = -j2 v - sm $- (33) 



whence one sees that 



-!- sin 6+ = -!- sin 0_. (34) 

J 4- J — 



This is precisely Snell's law of refraction, where we identify J_ 1 , J, 1 as the refractive indices of the materials 
on either side of the wall. Hence, total internal reflection occurs if 

|0_| > sin- 1 (J_/J+). (35) 

Note that l|34|l is independent of as it must be, since the sharp domain wall J{x) is invariant under 
dilation of the physical plane (x,y) <— ► (Xx,Xy). It is straightforward to solve the geodesic equations for 
g = f x (a\){da\ + da^) numerically and confirm the behaviour predicted. Figure El depicts bubble scattering 
off a sharp domain wall with J + /J_ = 2 at a range of incidence angles. 




Figure 2: Scattering of bubbles of width \ = 0.5 (solid curves) and width \ = 0.1 (dashed curves) off a domain 
wall, in the geodesic approximation. The region to the left of the vertical line x = has J = J_ = 1, while that to 
the right has J = J + = 2. Note that bubbles are reflected if their incidence angle exceeds 30°, in accordance with 
Snell's law. 

Although we have derived Snell's law within the framework of the geodesic approximation, one can rederive 
it using only field theoretic conservation laws. By Noether's theorem and t and y translation invariance, the 
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full field theory conserves E = T + V and 



p2 = I dxdy 7W^'^- (36) 

The key point is that the effective spacetime metric 77 is, except in a narrow strip containing the y-axis, the 
Minkowski metric, but with the "speed of light" differing on either side of the wall. Given a static solution 
strongly localized in one of the two flat regions, therefore, we can produce (almost) exact moving solutions by 
an appropriate Lorentz boost. Let no(x,y) be the static bubble given by u(z) — x The boosted bubble 
with velocity v- (cos 9- , sin 0_ ) , 

n_(t, x, y) = n (7- (x cos 6L — ysin6>_ — v_t), xsinO- + y cos6L), (37) 

where the contraction factor is 7_ = (1 — v^/J 2 .) - ^, tends to an exact solution as t — ► —00. Using this as 
initial data at t = —00, let us consider the field's subsequent evolution. Presumably the bubble impinges on 
the domain wall. Let us assume that the bubble is transmitted, and that its interaction with the wall is elastic, 
so no energy is lost to radiation. Then, as t — > 00 it is once again a free bubble moving at constant velocity, 
so must approach some translate of 

n_|_(t, x, y) = no (7+ (a; cos 9 + — y s'm9 + — v + t),x sin9 + + y cos 9 + ), (38) 

with contraction factor 7+ = (1 — v^/J 2 .) - ^. It is straightforward to compute E and P2 for n± as t — > ±00, 
yielding the conservation laws 

v + 7+ sin 9+ v_7_sin(9_ 
A/ 7+ = M 7 _ and + ,+ j2 -M = ] — 2 M, (39) 

where M = 4tt is again the rest energy of the bubble, which together imply (|34(l . Note that M is independent 
of Xj so this argument holds even if the interaction with the wall somehow changes the bubble's width. The 
limiting factor on the applicability of l|34|l is thus the assumption of elastic scattering off the wall, which one 
expects to hold at small incidence velocities, rather than the validity of the geodesic approximation itself. 



5 Circular lenses 

It is interesting to consider the interaction of bubbles with curved domain walls. The observation that bubbles 
undergo refraction, with J~ 1 identified as the refractive index, immediately gives a simple "geometric optics" 
model of their trajectories, which should hold whenever the bubble width is large relative to the domain wall 
width, but small relative to its radius of curvature. This raises the interesting possibility that lenses can be 
constructed which could be used to focus bubbles onto a target area. Such a lens would consist of a thin region 
bounded by a pair of circular arcs of large radius, within which J is suppressed relative to its value outside the 
region. This is a rather complicated choice of J(x, y) for the purposes of computing g, so we shall here test the 
feasibility of antifcrromagnetic lenses more indirectly, by studying the scattering of bubbles off complete disks 
of suppressed J, for which g can again be computed explicitly. Clearly, such circular "lenses" are far from the 
thin lenses used for focussing in optics. Nevertheless, by comparing bubble trajectories incident on a circular 
lens with their geometric optics analogues, we can provide strong evidence that bubble focussing is possible. 

We shall again make a thin domain wall approximation, so that J{z) is piecewise constant, taking one 
value on some disk and another value elsewhere. By conformal invariance we can, without loss of generality, 
assume the disk has unit radius. Hence 

^ - { t !?i > : < 4 °» 

with J_ < J + . We can again use formula l)28[l. but now 5_ C S 2 is the region on which \\u + a\ < 1. This is 
again a disk, whose stereographic image is bounded by the circle of radius x _1 centred at —a/x- By rotating 
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around the axis, this disk can be centred on the positive 112 axis, so that the segment of this axis from 
u 2 = (M — l)/x to U2 = (|a| + l)/x becomes a great circular arc in 5_ intersecting dS- orthogonally. Hence 



diam(iS_) = 



(kl+i)/x 



dt 



which, together with (|3f)|> and (|28J) gives 



l + < 2 



/x( a ) = 



1 - 



X 



1 



= 2 (tan" 



2tt 



1 



— tan 



F 



Vx 4 + 2x 2 (H 2 + l) + (|a| 2 -l) 5 



(41) 



(42) 



Note that f x depends on \a\ alone, is smooth, and has lim^i^^ f x (a) = AttJ^ 2 , as one would expect. The 
dependence on the bubble width x cannot, as in the case of a straight sharp domain wall, be scaled away. Unless 
X *C 1, therefore, we should expect bubble trajectories to diverge significantly from there optical analogues. 

For x sufficiently large x > Xi (~ 0-17 if J+/J-—2), (M*,g) may be isometrically embedded as a surface of 
revolution in R 3 , using the method described in 0. Figure shows numerically generated generating curves 
for these surfaces of revolution for various values of %: the surfaces themselves are swept out as the curves 
are rotated about an axis vertical in the page. It is not surprising that the isometric embedding fails for x 
small. Recall that in the limit % — > 0, g tends to a multiple of the metric on the physical plane, which for 
this choice of J is clearly singular. It can be roughly visualized as the flat plane with the unit disk cut out, 
expanded by a linear factor of J+/J_, its edge then re-attached to the hole. The surfaces depicted in figure EH 
are on the way to this geometry as x gets small. An interesting point to note is that a pinched neck develops 
in the surface for x < X2 (~ 0.36 for J + /J_ = 2). The two circles swept out by the points where the tangent 
to the curve is vertical are both geodesies. Hence there is a bifurcation in the geodesic flow at x = X2, & pair 
of periodic geodesies appearing as x is reduced through X2- Although we have deduced the existence of these 
periodic geodesies visually, it is easy to show that they exist for all x < X2 (even if x < Xi 1 so that no isometric 
embedding exists) by appealing to conservation of energy e and angular momentum pg = / x (|a|)|et| 2 #, where 
a = |a|e l61 . The corresponding bubble trajectories are circular orbits centred on z = 0, one inside the disk and 
one just outside. Such trajectories have no optical analogue. 




z=0.2 



Figure 3: Generating curves for isometrically embedded as a surface of revolution in M 3 for various x< ,n t ne 
case of a sharp disk inhomegeneity with J_ = 1, J + = 2. The surfaces are swept out as these curves are rotated 
about the axis x = 0. Note that a neck develops as x becomes small. 

It is again straightforward to solve the geodesic problem numerically. Figure 01 shows the scattering of 
bubbles initially moving along parallel trajectories for two different choices of x- A focussing effect is apparent, 
and the trajectories look optical in the case of very narrow bubbles. To make a quantitive comparison, we 
present in figure^lplots of the exit angle of the bubble as a function of its impact parameter b = linii^^oo ai (t) 
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for a variety of widths Xi i n comparison with the exit angle for light rays, 



7 ii g ht(&) = 2 ( sin 1 -j-b - ■ 



(43) 



Clearly, bubble scattering is increasingly optical as they become narrow, as one expects. Even for moderately 
narrow bubbles the geometric optics picture works well except for very glancing impacts (b close to 1). One 
area where the geometric optics model fails badly is scattering with b just exceeding 1. Light rays in this case 
miss the disk completely, so #ii g ht(&) = 0. Bubbles, on the other hand, can become almost trapped by the neck 
in M*, and scatter with nontrivial 9, having wound several times around the disk. The dependence of 9 on b 
is very sensitive in this regime, see figure 0] 

(b) 






J =1,J =2 



Figure 4: Scattering of bubbles off a circular "lens" with J + /J_ = 2: (a) and (b) show the expected focussing 
effect for \ — 0-3 and 0.01 respectively, while (c) shows the sensitive dependence of the trajectory on its initial data 
when it winds around the neck of (x — 0-3 again). 









J =1, J =2 

- + 






solid: % = 1, 1/10, 1/100 






dashed: light rays 
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Figure 5: The exit angle of a bubble trajectory scattering off a circular lens with J + /J_ = 2, as a function of 
impact parameter for x = 1. 0.1, 0.01 (solid curves, top to bottom), in comparison with the exit angle for light rays 
(dashed). 



If we choose J_ > J+, the disk causes bubble trajectories to diverge rather than focus, as one would expect. 
Again the trajectories are increasingly optical as x becomes small. A significant difference is that M \ can never 
be isometrically embedded in R 3 , and there are no periodic geodesies. 
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(a) (b) 




Figure 6: A circular bubble guide, analogous to a fibre optic cable. The region between the dashed circles has J 
suppressed by a factor of 2. Total internal reflexion keeps narrow bubbles (\ = 0.1) between the walls, (a). Even 
a bubble as wide as the track itself (x = 1) is guided along it, (b). 



6 Bubble guides 

The total internal reflexion phenomenon might be used to construct the analogue for bubbles of fibre optic 
cables: doped tracks of low J in a high J substrate which would guide bubbles along a prescribed path. A 
simple example, for which f x [a) can still be obtained exactly, is an annular track with 

'^ Z ^ | J + elsewhere 

where J_ < J+ as before. Once again (|28J) holds, but now S- (or rather its stereographic image) is the annular 
region centred on —a/x bounded by radii r±/x and r^fx- The area of this region can be computed in similar 
fashion to the previous cases (by finding the diameters of its inner and outer disks) , yielding the formula 

4-7T 1 1 

fxio) = jW+Hj2 - j2)[KHx,n) - h(\a\, X ,r 2 )}, (45) 

X 2 + Q? — r 2 

h(a,x, r) = — ; (46) 
y/ X 4 + 2x 2 (a 2 + r 2 ) + (a 2 - r 2 ) 2 

Bubble trajectories within such an annular bubble guide with J + /J_ = 2 and n = 9, r2 = 10, obtained 
numerically within the geodesic approximation, are depicted in figure A narrow bubble {x — 0.1, compared 
with the track width of unity) initially moving within the annulus is trapped by the domain walls and guided 
around it. Even a bubble as wide as the track (x = 1) is guided along it, although the bubble centre now 
strays beyond the domain walls, and the trajectory looks far from optical. 



where 



7 Conclusion 

We have argued that, in the near continuum regime, a two-dimensional isotropic Heisenberg antiferromagnet 
with inhomogeneous doping should be well modelled by the relativistic 0(3) sigma model on a spatially 
inhomogeneous spacetime. By conformal invariance, such a model supports static Belavin-Polyakov lumps, 
which we call bubbles. The trajectories of these bubbles are predicted to behave like light rays propagating in 
an inhomogeneous medium, with J -1 , the inverse exchange integral, identified with the refractive index. In 
particular, bubbles incident on a straight domain wall are refracted in accordance with SnelPs law, and total 
internal reflexion can occur. This result was derived both within the framework of the geodesic approximation, 
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and by appeal to field theoretic conservation laws. Its validity is not contingent on the validity of the geodesic 
approximation per se, but rather on the assumptions that the interaction of bubble and domain wall is elastic, 
which should certainly be true for low impact speeds, and that the scattering does not cause the bubble to 
collapse immediately (it must survive long enough to travel far from the domain wall for our conservation law 
argument to apply). Bubble collapse is an intense focus of current research, and it is still unclear whether 
collapse is generic. 

In the present context, it would be perverse to try to settle this issue by lattice simulation of the sigma 
model, since this is itself merely an approximation to the more fundamental spin lattice Numerical 
simulations of (|TJ seeking the refraction phenomena noted here are the logical next step, but lie beyond the 
scope of this paper. Even if bubble collapse does turn out to be a serious problem, there is still hope that bubble 
refraction may occur in real-life systems because, just as in the case of ferromagnets, higher order interactions 
may contribute extra terms which stabilize the bubble [U]|- These terms destroy the conformal invariance of 
the sigma model, so there would be corrections to Snell's law due to the J dependence of the bubble's rest 
mass, but the basic refractive behaviour may persist. The main obstacle to seeing such refraction in a real 
antiferromagnet is that, since neighbouring spins almost anti-align, bubbles have almost zero magnetization, 
and are thus difficult to observe experimentally. 

To illustrate the possibility, it is instructive to consider the spatially inhomogeneous system studied recently 
by Piette et al [TT]. 

S = J dtdxdy Ifyn ■ d»n + ^[(fyn ■ <%n)(cPn ■ d"n) - (fyn ■ d"n) 2 ] - ^ ^ (1 - n§), (47) 

which may be loosely interpreted as an antiferromagnet with homogeneous and isotropic exchange integral 
(J = 1), but with inhomogeneous on-site easy-axis anisotropy (the potential coefficient fj,(x,y) depends on 
position). This system supports approximate bubble type solutions, localized in regions of constant /x, but 
their width is no longer a free parameter, and both their width and rest-mass now depend on position. Let 
us assume again that we have domain wall doping, so fi(x, y) = /i_ for x < 0, fi(x, y) = /i+ > /i_ for x > 0. 
Then we expect the rest mass of a bubble in the right half-plane, M + , to exceed that of a bubble in the left 
half-plane, M_ (Piette et al find M + /M_ w 1.029 for \i_ = 0.6, /x+ = 0.7, for example). What happens when 
a bubble moving with speed u_ at incident angle 0- hits the domain wall? By energy conservation, it must 
be reflected if 7_M_ < M+, independent of 0_. Assume V- exceeds this threshold, v 2 _ > 1 — M^_/M^_, and 
the bubble crosses the barrier, escaping with speed v+ at angle + . Then energy conservation M_7_ = M + 7 + 
fixes v+ as a function of t>_, and clearly w + < u_. Again, the y-component of linear momentum is conserved, 
so w_7_ SXD.9—M— = w+7+ sin# + Af + , which together with energy conservation imply 

v+ sin 8 + = u_ sin . (48) 

Since w_ > v + refraction does occur (the bubble bends away from the normal) , and the bubble must undergo 
total internal reflexion if exceeds the v --dependent Brewster angle sin _1 [w + (y—)/v—\. 

Turning to theoretical high energy physics, the 0(3) sigma model has long been a favoured model for those 
interested in soliton dynamics in geometrically exotic spacetimes |S| ^1 ^3 El Q5I DEI ■ Our observation raises 
the possibility that some exotic spacetimes may be artificially engineered in the laboratory. Furthermore, since 
Snell's law emerges from conformal invariance (the soliton's rest mass is independent of J) and very general 
conservation law arguments, soliton refraction should not be a special feature of this particular model. One 
would expect it to be a feature of Yang-Mills instanton dynamics in 4 + 1 dimensions, for example. 
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